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Abstract 


The  ability  to  measure  the  roughness  of  surfaces  will  be  importara  for  general  purpose  machine  vision 
systems,  in  which  a  robot  attempts  to  passively  gain  irrformation  about  its  environment.  Although  the 
traditional  study  of  profilometry  is  well  developed,  its  methods  are  unsuitable  for  this  endeavor,  since 
they  require  specialized  apparatuses  such  as  laser  interferometers  and  styluses.  We  propose  a  method 
for  the  visual  determination  of  surface  roughness  which  does  not  require  strong  environmental  con- 
straints  such  as  coherera  light  and  small  samples  of  the  surface  in  question.  We  determine  surface 
roughness  by  the  measurement  of  the  sharpruss  of  the  edges  of  an  image  r^cted  in  a  rough  surface. 
We  choose  step  edges  since  they  will  be  present  in  many  environments,  and  because  the  methods  of  edge 
detection  and  localization  are  mature.  For  these  initial  experiments,  we  assume  a  common  surface 
roughness  rturdel  of  simple  form,  which  we  believe  models  milled  or  rolled  surfaces.  We  show  that  if  we 
also  assume  a  known  shape  for  the  reflective  surface  and  a  known  position  for  the  viewer,  we  may  calcu¬ 
late  a  roughness  parameter.  For  our  surface  roughness  model,  we  show  that  the  resulting  intensity  dis¬ 
tribution  may  be  given  in  closed  form.  We  then  show  the  results  of  experiments  performed  with  planar 
surfaces  in  a  controlled  environment  and  give  the  calculated  roughness  parameter.  We  compare  these 
results  with  the  values  of  the  roughness  as  measured  with  a  stylus  profilometer,  and  find  that  the  results 
are  quite  different.  We  discuss  the  possibilities  that  Ms  difference  is  due  to  the  non-Gaussian  character 
of  the  surface  height  distribution,  and  the  possibility  that  the  difference  arises  from  the  difference  in 
measurement  procedures.  We  next  show  some  extensions  to  other  type  of  roughness  functions,  and  dis¬ 
cuss  some  new  methods  of  roughness  estimation  which  may  yield  information  about  a  broader  class  of 
roughness  models.  Finedly,  we  discuss  the  additional  issues  in  extending  and  applying  this  approach  to 
a  general  vision  system. 


1.  Introduction 


The  detennination  of  surface  roughness  is  important  in  many  areas,  from  quality  control  in  factories,  to  general  pur¬ 
pose  machine  vision  systems.  Although  there  exist  many  instruments  which  measure  roughness  accurately,  such  as 
laser  interferometric  devices  and  stylus  profilometers,  these  instruments  all  require  strong  constraints  on  die  environ¬ 
ment  and  the  samples,  for  example,  the  necessity  for  coherent  light  or  limits  on  the  sample  size.  These  constraints 
make  such  devices  impractical  for  a  general  purpose  vision  system,  in  which  a  robot  attempts  to  gain  information 
about  its  environment.  It  is  obvious  that  we  can  obtain  some  information  about  roughness  without  these  kinds  of  con¬ 
straints,  since  humans  can  often  make  crude  guesses  of  surface  roughness,  and  can  even  differentiate  surfaces  acced¬ 
ing  to  roughness,  solely  on  the  basis  of  visual  information.  For  example,  technicians  in  mills  sometimes  check  die 
roughness  of  sheet  metal  rolling  equipment  by  observing  the  reflection  of  their  fingers  in  die  rollers.  As  another 
example,  consider  Figure  1 .  This  is  a  picture  of  a  stainless  steel  ball  bearing  that  has  been  roughened  by  sand  blasting. 
In  the  bearing,  we  see  imperfect  images  of  the  objects  in  the  room,  including  rows  of  fluorescent  lights,  windows, 
pens,  and  even  the  photographer.  Most  people  will  be  able  to  make  a  guess  at  the  roughness  of  the  metal,  and  say  that 


Figure  1:  An  image  with  the  reflections  of  step  edges 


it  is  “smooth",  “a  little  rough",  and  so  on,  even  though  they  lack  sophisticated  instruments.  We  believe  that  they  do 
this  by  analyzing  the  blurriness  of  the  reflected  image  of  the  environment  in  the  surface.  We  propose  to  mimic  this 
ability  by  analyzing  reflected  step  edges;  step  edges  are  valuable  for  analysis  because  they  provide  considerable  con¬ 
straint  on  the  form  of  the  reflected  image,  are  present  in  almost  all  environments,  and  will  be  relatively  easy  to  find  in 
images,  since  the  field  of  edge  detection  is  well  developed.  We  will  first  present  a  taxonomy  of  the  commonly  used 
surface  roughness  models.  Then,  in  order  to  describe  reflected  step  edges,  we  must  choose  one  of  these  models,  so  we 
use  the  one  which  is  most  common  among  computer  vision  researchers,  and  which  models  surfaces  which  have  been 
rolled  or  milled.  We  believe  this  model  is  general  enough  to  yield  useful  results.  We  then  find  the  form  of  the  image  of 
the  reflected  step  edge  as  a  function  of  roughness,  and  show  the  results  of  our  experiments  with  rough  surfaces. 


2.  Roughness  and  reflection 


Let  us  begin  by  describing  some  of  the  rough  surface  models  used  by  machine  vision  researchers.  Figure  2  presents  a 
simple  taxonomy  of  these  models;  this  listing  is  by  no  means  complete,  but  clarifies  the  relationships  between  some 


of  the  commonly  used  models.  We  may  parametrize  a  surface  by  means  of  its  height  values  or  its  slopes,  and  the 
choice  will  determine  if  we  want  to  solve  for  the  reflected  light  intensity  by  means  of  the  methods  of  physical  optics 
or  of  geometrical  optics.  The  diagram  shows  that  there  is  some  overli^  between  the  two  formalisms,  because  some 
height-based  models  have  a  slope-based  interpretation.  Next,  we  may  consider  the  directionality  of  the  roughness, 
although  most  machine  vision  researchers  consider  only  the  limiting  cases  of  isotropic  and  one-dimensional  rough¬ 
ness.  It  is  also  common  to  assume  that  the  heights  or  slopes  follow  a  Gaussian  distribution.  Note  that  the  Gaussian 
height  models  have  a  slope  interpretation,  and  thus  also  fall  into  the  category  of  slope-based  models.  Not  all  Gaussian 
slope  models  have  a  height-based  interpretation,  however.  The  isotropic  and  one-dimensional  Gaussian  models  are 
the  most  commonly  used  classes  of  models.  The  isotropic  Gaussian  height  model  and  the  one-dimensional  Gaussian 
height  model  are  those  proposed  by  Beckmann  and  Spizzichino[l].  The  isotropic  Gaussian  slope  model  is  that  pro¬ 
posed  by  Torrance  and  Spanow[12],  and  is  probably  the  most  conunonly  used  model  in  computer  vision 
research[3,4, 5,7,9].  The  one-dimensional  Gaussian  slope  model  is  die  one-dimensional  analogue  of  this  model.  Note 
also  that  the  isotropic  Gaussian  height  model  and  the  isotropic  Gaussian  slope  model  are  equivalent  in  the  limiting 
case,  as  are  the  one-dimensional  Gaussian  height  model  and  the  one-dimensional  Gaussian  slope  model[l,81.  Thus, 
these  classes  overlap.  Other  models,  such  as  fractal  height  models,  appear  in  the  optics  literature[6],  but  are  uncom¬ 
mon  in  machine  vision  research,  probably  due  to  their  complex  nature.  These  models  do  not  appear  in  the  diagram, 
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since  their  relationship  to  other  models  is  uncertain,  although  they  would  belong  in  the  category  of  height-based  mod¬ 
els. 

We  will  quantify  the  scattering  from  rough  surfaces  with  a  slope-based  model,  in  which  light  reflects  from  small,  mir¬ 
ror-like  facets  which  comprise  the  surface;  we  will  use  the  methods  of  geometrical  optics  to  calculate  the  reflected 
radiance  distribution.  The^  slope-based  models  do  not  deal  with  the  individual  facets  on  the  surface,  but  rather  the 
statistical  properties  of  the  set  of  all  facets  on  the  surface. 

Let  us  now  derive  the  radiance  of  light  reflected  by  a  rough  faceted  surface.  This  derivation  parallels  that  given  by 
Torrance  and  SparTow[12].  Consider  a  small  planar  patch  of  tl»  surface  under  study,  of  size  dA  (see  Figure  3 ). 


Assume  an  infinitesimal  source  of  radiance  L,.  subtending  a  solid  angle  d(0.  with  direction  of  incidence,  in  spherical 
coordinates,  (0,..<pp.  (The  normal  to  the  mean  surface  is  in  direction  (0,0).)  The  reflected  radiance  is  dL^,  with  direc¬ 
tion  (8^,<p p.  For  a  given  incidence  direction  and  reflectance  direction  only  those  facets  whose  normals  bisect  the 
angle  between  these  directions  will  be  able  to  reflect  light  from  the  source  to  the  receiver.  Let  the  direction  of  the  nor¬ 
mals  to  these  facets  be  (Oytpp.  Let  us  denote  the  local  angle  of  incidence,  or  the  angle  between  the  source  direction 

and  the  direction  of  the  normals  to  the  privileged  facets,  as  6.^. 

Now,  in  order  to  calculate  the  radiance  of  the  patch,  we  must  know  the  total  area  of  the  privileged  facets,  which  will 
reflect  light  from  the  source  to  the  receiver.  We  therefore  define  a  function,  S  (9^,  tpp  ,  the  facet  area  (frstribution  func¬ 
tion  (FADF),  which  represents  the  total  area  of  the  privileged  facets  per  unit  solid  angle  subtended  by  the  normals  to 
these  facets  per  unit  area  of  the  mean  surface.  Thus,  the  area  of  the  facets  present  on  an  infinitesimal  mean  surface 
patch  of  area  dA  with  normals  in  the  infinitesimal  solid  angle  day  is  5(6^  <pp  dfo^dA.  The  power  incident  on  the 
privileged  facets  is  then 

/>.  =  Z,.</©.cos0.^5(0^<ppdoydA  (1) 

This  light  will  then  reflect  off  the  facets  and  into  the  receiver,  and  may  be  attenuated  upon  reflection.  From  the  point 


of  view  of  the  receiver,  the  power  received  from  the  facets  may  be  written  as 


=  FGP^  =  dL^cosQ^dadA  (2) 

where  dL^  is  the  radiance  of  the  patch  as  seen  from  the  viewing  direction,  and  F  is  the  Fresnel  coefficient  for  the  sur¬ 
face.  The  Fresnel  coefficient  will  be  a  function  of  the  local  angle  of  incidence,  8^.^.  Some  researchers[]  also  model  the 

attenuation  due  to  the  shadowing  of  one  facet  by  another,  and  quantify  this  with  G,  the  geometrical  attenuation  factor. 
It  is  a  function  of  the  incident  and  exitant  directions.  Torrance  and  Sparrow  provide  a  derivation  of  this  factOT.  Substi¬ 
tuting  the  expression  for  the  incident  power  into  the  equation  for  the  reflected  power,  we  see  that  the  radiance  of  the 
surface  patch  is 


dL 


r 


L.</(o.cos8. ^FGS(0^  qy) 
cosFd©^ 


Most  researchers  use  the  simplifying  relation[8,12 


da 

da,  =  - 

/  4cos6j.^ 

We,  however,  will  most  often  use  the  relation 


(3) 


(4) 


da.  =  da^  (5) 

We  have  found  that  numerical  methods  of  calculating  the  reflected  radiance  often  perform  better  if  we  use  this  rela¬ 
tion.  Additionally,  the  equation  for  the  reflected  radiance  has  the  same  form  for  both  the  one-dimensional  and  two- 
dimensional  cases  if  we  make  this  substitution. 


For  an  extended  source,  the  reflected  radiance  is  then 


Lj.FGcos0^.^5  (0^  qy)  da^ 


(6) 


where  the  integral  is  over  the  hemisphere  for  which  0^  is  less  than  - . 

There  is  also  a  normalization  condition  which  any  realistic  FADF  must  satisfy.  Consider  one  facet  of  the  surface  and 
let  its  area  be  dF.  The  area  of  the  projection  of  this  element  onto  the  mean  surface  is 

dFcosQy  (7) 

If  we  were  to  add  up  all  the  projected  areas  of  all  facets  on  the  surface,  we  would  obtain  dA ,  the  area  of  the  patch,  i.e. 

dA  =  JdFcos0^ 

We  also  know  that 


dF  =  S(8^qy)d©^A 


so  that 


dA  -  Js  (0^  qy)  cosQ^a^A 
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and 

;s(e^V  cosGydcy  =  1 

which  is  the  normalization  condition. 

3.  One-dimensional  Gaussian  heights 


(8) 


3.1.  Theory 

The  results  so  far  pertain  to  any  slope-based  model,  or  FADE  In  (wder  to  proceed  further,  we  must  now  choose  a  spe¬ 
cific  model  for  our  rough  surface.  Of  particular  interest  is  oi»  which  we  call  the  “one-dimensional  Gaussian  heighf  ’ 
distribution  function.  We  believe  this  function  approximately  models  surfaces  which  have  been  rolled  or  milled.  It  is 
also  one  of  the  more  commonly  used  models.  Consider  a  surface  which  varies  in  height  along  one  direction,  but 
which  is  constant  in  height  along  the  perpendicular  direction.  The  surface  thus  consists  of  a  series  of  parallel  grooves. 
Assume  that  the  height  values  along  the  first  axis  are  random  and  follow  a  Gaussian  distribution.  We  wish  to  model 
this  surface  as  a  collection  of  facets.  As  shown  by  Beckmann  and  Spizzichino[l],  such  a  surface  is  described  by 


5(0p  = 


(9) 


where  P  is  the  roughness  parameter,  and  6^  is  the  angle  between  the  normal  to  the  surface  patch  and  the  facet  normal. 


All  the  facet  normals  lie  in  a  plane,  and  the  standard  deviation  of  the  facet  slopes  is 


B 

-= .  Note  that  for  this  case,  doo. 


is  replaced  by  dQy  in  equation  (6),  since  all  the  facets  are  coplanar  and  may  be  described  by  one  variable.  As  noted  by 


Beckmann  and  Spizzichino[l]  and  Nayar[8],  the  function  may  be  approximated  by  a  Gaussian  in  the  limit  of  sntall  P, 
and  thus  the  model  of  Torrance  and  Sparrow,  with  its  Gaussian  distribution  function,  is  a  reasonable  approximation  to 
this  function  for  the  case  of  a  fairly  smooth  surface.  This  facet  area  distribution  function  is  more  realistic  than  the 
Gaussian  model,  however,  since  it  satisfies  the  normalization  condition  (equation  (8)). 


Now,  let  us  determine  the  radiance  pattern  for  the  reflection  of  a  step  edge  in  a  rough  surface  with  the  given  facet  area 
distribution  function.  We  will  perform  our  calculations  in  the  coordinate  system  of  the  surface  patch.  Therefore,  let 
the  normal  to  the  patch  be  along  the  z-axis,  as  shown  in  Figure  2,  and  let  the  grooves  lie  parallel  to  the  y-axis  (i.e.  the 
height  is  constant  for  constant  x).  Let  the  source  be  a  Lambertian  half-plane  of  radiance  /.  Although  this  may  seem  to 
be  a  very  restrictive  assumption,  because  we  are  studying  the  reflected  images  in  the  regions  local  to  edges,  it  is  not 
excessively  restrictive;  if  we  wish  to  study  the  interaction  of  neighboring  edges,  we  will  have  to  drop  the  half-plane 
assumption.  Let  us  first  consider  the  easiest  case,  in  which  the  view  vector  lies  in  the  plane  of  the  surface  normals, 
i.e.,  the  x-z  plane.  Assume  the  edge  of  the  source  is  parallel  to  the  y-axis.  We  will  remove  these  restrictions  shortly. 
As  we  see  in  Figure  4,  light  is  incident  on  the  patch  through  a  range  of  values  of  the  angle  6^. .  Let  0^.^  denote  the  start¬ 
ing  value  and  0 .  j  the  ending  value  of  this  range.  We  also  know  the  direction  of  the  view  vector,  and  hence  0^.  for  the 
patch.  Because 


0.  +  0^ 
2 


(10) 


we  therefore  know  the  range  of  values  of  0^  for  which  light  will  reflect  from  the  source  into  the  camera.  Let  us  denote 


the  starting  and  ending  values  of  this  range  by  0^  and  9^ ,  respectively.  Now,  let  us  return  to  equation  (2),  which 

gives  the  reflected  radiance  for  a  patch.  We  will  make  the  simplifying  assumptions  that  F  =  1  and  C  =  1 .  The  first 
assumption  bolds  fairly  well  over  a  iaige  range  of  angles  for  most  metals.  The  second  assumption  holds  when  both  6. 

and  6^  are  less  than  approximately  43  degFees[8].  We  expect  some  error  due  to  this  assumption,  but  not  a  substantial 
amount  With  these  assumptions,  the  reflected  radiance  will  be  given  by 


rtan0,\^ 


^^/cos9.,e 

‘■r-l 


-T 


d0. 


9^  p.,^(cOS0y.)^COS0^ 

Because  the  angle  at  which  light  is  incident  on  a  facet  must  equal  the  angle  at  which  light  exits  the  facet 
cosO.j^  =  cos  (9.- op  =  cos  (0^-0^)  =  cos0^os0^  + sin0^in0^ 

and  thus 


(11) 


I 


'’/o 


,~vT ) 


1800,^2 


le  ^  '  dO.  le 

^  +  — 


tan0^tan0^d0y. 


P Tit  (cos 9^)'' 


p7jt(cos0p‘ 


(12) 


and 
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The  radiance  of  the  reflected  image  will  change  from  point  to  point,  since  the  values  of  0^,  0^ ,  and  0^  will  change 
from  point  to  point  on  the  surface. 

Now  let  us  remove  the  assumptions  that  the  edge  of  the  source  is  parallel  to  the  grooves  of  the  rough  reflecting  sur¬ 
face,  and  that  the  view  vector  lies  in  the  x-z  plane.  We  note  that  the  rough  surface  is  still  parametrize  by  a  single  vari¬ 
able,  0^  so  that  we  may  still  calculate  the  reflected  radiance  by  means  of  equation  (11).  Note,  however,  that  for  this 

case,  cos0.^  is  given  by 


cos0.^  =  cos0^os0^+ sin0^in0^cos((p^-(p^) 
Because  the  surface  is  “one-dimensionally“  rough,  however,  qy  =  0,  and  thus 

cos0.^  =  cos0^os0^-e  sin0^in0^cos<p^ 


so  that 


(14) 


B 

—  tan0^cos<p^ 


^  ^  -e 

1  p  J 

V 

(15) 


This  description  of  the  reflection  of  light  from  rough  surfaces  is  not  complete,  however.  As  shown  by  Beckmann  and 
Spizzichino,  there  is  a  purely  specular  component  of  reflection  which  can  only  be  explained  by  considering  the  wave 
nature  of  the  light.  Nayar[8]  calls  this  component  the  “specular  spike  component  of  reflection”,  and  that  described  by 
equation  (15)  as  the  “specular  lobe  component”.  Apart  from  a  multiplicative  constant,  the  spike  component  of  reflec¬ 
tion  is  the  same  as  the  reflection  that  would  occur  if  the  surface  were  completely  smooth.  Thus,  we  model  the  light 
reflected  from  the  rough  surface  as 


L,  =  A(^)-hB 
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where  A  and  B  describe  the  relative  intensities  of  the  two  components,  and  (^)  means  that  the  intensity  is  non-zero 
only  if  the  light  may  reflect  specularly  from  the  source  to  the  receiver. 


3.2.  Experiment 

We  tested  the  correctness  of  the  theory  by  calculating  the  roughness  parameter  for  several  surfaces.  The  samples  were 
circular  stainless  steel  disks,  three  inches  in  diameter,  which  had  been  milled  to  produced  a  one-dimensionally  rough 
surface  on  one  face  of  each  disk.  The  values  of  the  arithmetic  average  roughness,  [11],  as  stated  by  the  manufac¬ 
turer  were  1.5,  3.0, 4.5,  7.5,  and  9.5  microinches.  The  camera  was  a  12-bit  Photometries  Star  1  camera.  A  large  piece 
of  white  paper  was  used  as  the  step  edge  source.  Because  the  source  had  a  much  greater  area  than  the  steel  plates,  it 
was  an  appropriate  choice  to  approximate  a  step  edge.  Figure  5  shows  the  experimental  conflguration.  The  distance 
from  the  disks  to  the  camera  was  approximately  120cm,  and  the  distance  from  the  plane  of  the  samples  to  the  plane  of 
the  source  was  9.57cm.  The  viewing  angle,  0^,  was  7.29  degrees. 


r 


t(q)  view 


front  left  view 


plate 


camera 


Figure  5:  Schematic  of  the  experimental  configuration 


In  order  to  apply  the  results  of  the  theory,  the  intensity  of  the  reflected  step  edge  should  be  studied  along  the  direction 
perpendicular  to  the  reflected  edge.  For  a  generic  configuration  of  source,  reflective  surface,  and  camera,  it  would 
therefore  be  necessary  to  find  the  reflected  edge  and  determine  its  direction.  This  is  not  difficult  in  general,  and  would 
be  the  first  step  in  a  roughness  estimation  vision  system.  We  felt  that  the  extra  calculation  involved  in  edge  detection 
would  add  nothing  to  the  test  of  the  theory,  and  thus  the  reflected  edge  was  aligned  vertically  in  the  image.  Similarly, 
the  grooves  of  the  rough  surface  lay  along  this  direction.  This  was  done  merely  to  simplify  the  calculations,  and  in  no 
way  represents  a  loss  of  generality  in  the  test  of  the  theory.  Note  that  since  we  assume  an  orthographic  projection, 

is  a  constant  for  all  points  on  the  reflective  surface.  For  this  configuration,  <p^  -  0,  and  equation  (IS)  reduces  to  equa¬ 
tion  (13).  We  may  derive  all  the  necessary  terms  in  this  equation  from  Figure  6.  First,  let  us  find  the  values  of  6^  and 
0^  for  every  point  on  the  steel  plate.  Let  be  a  point  on  the  source,  and  be  a  point  on  the  object  We  see  that 


Q.  =  atan  ( 


D  ’ 


whCTC  D  is  the  distance  between  the  source  plane  and  the  object  plane.  Also, 


®/  = 


e,.+0^ 


atan  ( 


^src  ^obj 


)+e^ 


We  therefore  know  the  initial  and  final  values  of  0.  for  all  points  on  the  object  and  we  may  therefore  use  equation 
(13). 


The  images  used  as  data  appear  in  Figure  7  to  Figure  9  in  order  of  increasing  roughness.  The  appearance  of  the 
reflected  step  edges  is  as  one  would  expect  from  intuition. 
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Figure  6:  The  reflection  geometry  for  the  experimental  conflguration 


Figure  7:  The  reflected  image  of  the  step  edge  in  the  microinches  sample  and  the  3.0  microinches 

Rj^  sample 


Now,  in  order  to  fit  the  model  to  the  data,  we  must  account  for  several  effects  which  change  the  appearance  of  the 
reflected  image,  but  which  are  not  included  in  equation  (13).  First,  note  that  the  dimmest  parts  of  the  image  do  not 
have  a  zero  intensify,  but  rather  there  appears  to  be  a  constant  offset  due  to  ambient  light  and  perhaps  camera  effects. 
We  model  this  with  the  constant,  C.  Next,  note  that  the  value  given  for  0^  is  the  result  of  a  measurement  and  is  there- 


Figure  8:  The  reflected  image  of  the  step  edge  in  the  4^  microinches  sample  and  the  IS  microinches 

sample 


Figure  9:  The  reflected  image  of  the  step  edge  in  the  9.5  microinches  R^  sample 
fore  not  exact;  the  same  is  true  for  0..  Apart  from  the  term  tan 9^,  which  is  a  small  number  for  this  experimental  con- 
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figuration  and  which  is  multiplifd  by  other  small  numbers,  0^  and  0.  only  tq>pear  as  part  of  their  sum.  Therefore,  we 
may  model  the  errors  in  both  0^  and  0.  by  replacing  0^  with  0^  +  A  0  .  The  final  model  for  the  reflected  edge  is  then 


r  1 

e 

1 

1 

! 

:ca 

1  p  J 

< 

)) 

+  C  (17) 


where  the  first  component  was  Al  for  0^.^  <  0^  <  0^.  j ,  and  0  otherwise.  This  term  models  the  specular  spike  compo¬ 
nent  of  reflection  and  is  simply  the  reflected  image  which  would  have  occurred  had  the  surface  been  mirror-like.  The 
constant  B  gives  the  strength  of  the  specular  lobe  component,  and  the  other  terms  are  as  described  above.  We  allow 
the  constant  /  to  equal  one,  since  the  terms  A  and  B  can  model  any  changes  in  the  value  of  /.  Our  goal  is  to  find  the 
best  fit  values  of  A ,  B,  C,  p,  and  A  0  for  the  data. 


In  order  to  mitigate  the  effects  of  noise,  the  scanlines  of  ettth  image  were  averaged  to  yield  a  representative  reflected 


3.0  microinches  roughness 

Figure  11:  The  average  of  tbe  scanlines 


the  parameters  were  found  as  follows.  First,  the  minimum  value  of  the  reflected  edge  was  taken  as  the  initial  guess  for 
C.  Next,  a  step  edge  with  minimum  value  C  and  maximum  value  equal  to  the  maximum  value  of  the  reflected  edge 
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Figure  12:  The  average  of  the  scanlines 
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Figure  13:  The  average  of  the  scaniines 

was  formed.  The  sum  of  squared  differences  between  this  edge  and  the  reflected  edge  curve  was  found  for  different 
values  of  A  9.  The  minimum  of  the  sum  of  squared  differences  with  respect  to  A  6  was  found  by  bracketing  and 
bisection,  and  the  value  of  A  6  which  gave  this  minimum  was  taken  as  the  best  fit  value  of  A  9.  Next,  the  pixel  at 
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Flgnre  14:  Tlw  avmite  of  thetcaidiBn 

which  the  discontinuity  in  the  step  edge  occuned  was  found.  The  intensities  a  few  (usually  three)  pixels  on  eitho'  side 
of  the  central  pixel  woe  recorded,  and  their  difference  was  taken  as  the  initial  value  of  A .  The  first  guess  for  B  was 
then  the  maximum  value  of  the  reflected  edge,  minus  C,  minus  A . 

Once  the  initial  values  far  A,  B,  and  C  were  determined,  the  sum  of  squared  differences  between  die  predicted  and 
measured  curves  could  be  found  for  a  given  p,  and  the  best  fit  value  of  p  could  be  found  by  gradient  descent  When 
this  reached  a  minimum,  singular  value  decomposition  was  used  to  find  die  best  fit  values  of  A,  B,  and  C  for  this 
value  of  The  algoridun  alternated  between  fitting  for  P  and  fitting  for  A,  £,  and  C  until  convergence  was  reached. 
The  best  fit  values  for  each  sample  are  shown  in  Table  1;  as  we  noted  eariier,  we  may  calculate  die  standard  deviation 
for  the  facet  area  distribution  foncdon  from  3,  and  the  last  column  of  Table  1  shows  these  values  in  units  of  degrees. 
The  corresponding  curves  and  the  data  are  shown  in  Figure  IS  through  Figure  19.  The  solid  lines  represent  the  fitted 
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curve,  and  the  dashed  lines  show  die  average  scanline  plus  and  minus  the  error  bars. 

Best  fit  values  of  the  parameters 


sample 

A 

B 

C 

AO 

P 

std.  dev 

1.5  microinches 

2300 

1290 

171 

-0.0272 

0.0770 

3.12 

3.0  microinches 

1220 

2670 

251 

-0.0312 

0.0850 

3.44 

4.5  microinches 

298 

3510 

373 

-0.0440 

0.102 

4.13 

7.5  microinches 

311 

3360 

299 

-0.0368 

0.0990 

4.01 

9.5  microinches 

109 

3750 

318 

-0.104 

0.210 

8.51 

Ihbiel 


1.5  microinches  roughness 

Figure  15:  The  best  fit  scanline  and  error  bars 

We  see  that  the  fitting  process  orders  the  samples  by  roughness,  except  for  the  7.5  microinches  roughness  sample.  We 
will  discuss  possible  reasons  for  this  error  shortly.  The  values  found  for  P  are  reasonable,  and  are  similar  to  those 


4.5  microinch^  roughness 

Figure  17:  The  best  fit  scanline  and  error  bars 

found  by  other  computer  vision  researchers[9].  Although  the  fitting  process  gives  results  consistent  with  our  intuition, 
we  must  check  the  goodness  of  the  fit,  and  compare  the  results  with  roughness  measurements  made  by  other  methods. 


7.5  mlerolnchea  roughness 

Figure  18:  The  best  fit  scanline  and  error  tars 


9.5  microinches  roughness 

Figure  19:  The  best  fit  scanline  and  error  bars 


The  best  fit  curves  were  compared  to  the  data,  and  the  statistics  quantifying  the  goodness  of  fit  are  shown  in  Table  2. 


Chi  squared  value  for  the  fit 


sample 

r 

1.5  microinches 

174191 

0.000000 

3.0  micFoinches 

398126 

0.000000 

4.5  microinches 

181983 

0.000000 

7.5  microinches 

222623 

0.000000 

9.5  microinches 

87473 

0.000000 

Table! 

Although  the  curves  appear  reasonable,  but  not  exact,  the  fits  are  not  acceptable.  Note  that  the  oror  bars  are  probably 
too  latge,  since  we  averaged  over  many  scanlines,  and  thus  included  variations  in  the  source  intensity  in  our  calcula¬ 
tions  of  the  error.  Thus,  the  fits  are  certainly  unacceptable  for  the  more  realistic  case  in  which  the  errors  are  smaller. 
There  are  several  possible  reasons  for  the  poor  quality  of  the  fits.  First,  the  model  of  the  reflected  light,  equation  (13), 
is  technically  incorrect  Equation  (13)  should  only  be  used  for  very  rough  surfaces  (those  with  P  »  1 )  [1].  In  this  case, 
there  should  be  no  specular  spike  component  and  thus,  for  all  our  samples  which  showed  an  appreciable  spike  com¬ 
ponent  equation  (13)  is  not  tqrpropriate.  Notice  that  the  value  of  is  lowest  for  the  9.S  microincbes  sample,  pre¬ 
sumably  because  the  model  is  more  appropriate  for  the  rougher  surface.  Another  possible  source  of  error  is  the 
chance  that  the  surface  has  a  non-Gaussian  distribution  of  heights.  We  will  discuss  this  possibility  after  we  compare 
our  results  with  roughness  measurements  made  by  other  methods. 

Although  our  method  produces  roughness  estimates  comparable  to  those  of  the  methods  of  other  computer  vision 
researchers,  we  wish  to  determine  if  the  calculated  roughness  corresponds  to  the  roughness  for  the  surface  as  mea¬ 
sured  by  conventional  means.  Before  we  describe  our  experiments,  however,  we  mention  a  caveat  As  noted  by  Ben¬ 
nett  and  Mattson[2],  different  values  of  surface  roughness  will  be  obtained  for  the  same  surface  with  different 
instruments.  As  a  result  stylus  profilometers  generally  measure  higher  values  for  the  roughness  of  a  surface  than  do 
optical  profilometers.  Keeping  this  in  mind,  we  performed  the  following  measurements.  A  Tencor  Alpha  Step  200 
stylus  profilometer  was  used  to  measure  the  surface  profiles  of  the  samples.  For  each  scan,  the  profilometer  measured 
the  surface  height  at  2000  points  over  a  400  micrometer  distance,  and  ten  scans  were  made  for  each  sample.  The 
approximation  of  the  surface  by  facets  is  a  difficult  problem;  however,  for  a  surface  with  a  Gaussian  distribution  of 
heights,  we  know  that  the  facet  area  distribution  function  for  such  a  surface  is[l] 

r^(tanep^ 

5(6)  =  - ^ - -e 

^  laJnicosQy)^ 


(18) 
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where  CT  is  the  rms  height  and  T  the  correlation  length  of  the  surface.  Therefore,  ^  ■  The  calculated  nns 

roughness  values  appear  in  Table  3.  The  values  in  microinches  are  given  for  comparison  with  the  stated  roughness 


Measured  rms  roughness  value 


sample 

rms 

micrometers 

roughness 

microinches 

1 . 5  microinches 

0.0366  ±0.0002 

1.44 

3 . 0  microinches 

0.0690  ±0.0003 

2.72 

4 . 5  microinches 

0.136  ±0.001 

5.34 

7 . 5  microinches 

0.241  ±0.001 

9.48 

9 . 5  microinches 

0.255  ±0.001 

10.05 

IhbleS 


values.  Remember  that  the  roughness  values  stated  by  the  manufacturer  wne  the  arithmetic  roughness,  rather  than  the 
root  mean  square  roughness,  and  therefore  need  be  the  same  as  the  values  we  calculate.  Notice  the  unusually  high 
value  for  the  7.5  microinches  roughness  sanqtle.  Note  that,  except  for  the  fourth  sample,  the  nos  values  place  the 
samples  in  the  same  order  as  did  the  edge  fitting  method,  and  simple  intuition.  Next,  the  correlation  lengths  of  the  sur¬ 
faces  were  found  by  taking  the  autoconelation  of  each  profile,  and  then  fitting  a  function  of  the  fonn 


in  the  least  squares  sense  to  each  of  the  autocorrelation  functions.  The  correlation  lengths  for  the  surfaces  tq>pear  in 
Table  4.  Again,  notice  the  large  value  for  the  7.5  microinches  roughness  sample.  Graphs  of  the  scans  from  die  profilo- 
meter  confirm  that  this  sample  has  a  much  lower  fundamental  firequency  for  its  siuface  than  do  the  other  samples. 
These  two  high  values  could  explain  why  the  facet  area  distribution  function  for  this  sample  is  almost  the  same  as  that 
for  the  4.5  microinches  sample,  rather  than  showing  a  broader  distribution.  The  standard  deviations  of  the  slope  dis- 


Calculated  correlation  length 


sample  correlation  length  (micrometers) 


1 . 5  microinches 

5.8  ±2.7 

3 . 0  microinches 

8.6  ±3.8 

4 . 5  microinches 

5.1  ±1.2 

7 . 5  microinches 

21.2  ±2.4 

9 . 5  microinches 

4.2  ±0.8 

Table  4 

tributions  are  shown  in  Table  5.  As  can  easily  be  seen,  the  values  found  by  this  method  do  not  correspond  at  all  to 

Standard  deviation  of  the  slope  distribution 

sample 

standard  deviation  (degrees) 

1.5  microinches 

0.509  ±0.237 

3 . 0  microinches 

0.652  ±0.230 

4 . 5  microinches 

2.14  ±0.61 

7 . 5  microinches 

0.920  ±0.149 

9 . 5  microinches 

4.86  ±  1.22 

Table  5 


those  found  by  the  optical  method  and  listed  in  Table  3.  This  discrepancy  could  be  due  to  the  poor  quality  of  the  esti¬ 
mates  of  the  correlation  length,  or  to  the  fact  that  two  different  methods  of  measurement  were  used  for  these  calcula¬ 
tions.  Also,  it  is  possible  that  the  surfaces  do  not  exhibit  a  Gaussian  distribution  of  heights,  and  that  the  facet  model  of 
reflection  is  not  tqiplicable  to  these  samples.  We  study  this  possibility  next. 


page  21 


Figure  20  to  Figure  24  show  the  height  histograms  for  the  five  samples,  as  determined  by  the  prohlometer  data. 


micmmeters 


1.5  micFoinches  sample 

Figure  20:  Number  of  occurrences  of  each  height  value 


3.0  tnicroinches  sample 

Figure  21:  Number  of  occurrences  of  each  height  value 

Except  for  the  7.5  microinches  sample,  the  distributions  appear  reasonably  Gaussian.  We  must,  however,  test  to  see  if 
this  guess  is  correct.  The  number  of  measured  height  values  that  fall  within  a  certain  range  is  given  by  a  binomial  dis- 


-1.500  -1.000  -0.500  0.000  0.500  1.000  1.500 

micrometers 


4.5  microinches  saiiq>le 

Figore  22:  Number  of  occurrences  of  each  height  value 
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micrometers 

7.5  microinches  sample 

Figure  23:  Number  of  occurrences  of  each  height  value 


tribution[10].  We  therefore  model  these  histograms  as  binomial  distributions,  and  find  the  values  of  the  mean  height 
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micrometers 


9.5  microinches  sample 

Figure  24:  Number  of  occurrences  of  each  height  value 

and  roughness  which  give  the  best  lit  to  the  data.  These  values  appear  in  Table  6.  Notice  that  the  mean  values  arc  not 

Mean  height  and  roughness  parameter 


sample 

mean  height  (micrometers) 

P 

1.5  microinches 

-0.0025 

0.0588 

3.0  microinches 

-0.0025 

0.104 

4.5  microinches 

-0.0075 

0.216 

7.5  microinches 

-0.0025 

0.351 

9.5  microinches 

-0.0175 

0.443 

Table  6 


zero,  as  they  should  be,  and  that  the  values  of  the  roughness  parameter  are  not  the  same  as  those  shown  in  Table  3. 
This  is  an  indication  that  the  surfaces  do  not  follow  a  Gaussian  height  distribution.  We  next  perform  chi  squared  test¬ 
ing  to  determine  if  the  fit  is  justified.  The  results  of  this  procedure  show  that  the  surfaces  actually  are  not  Gaussian. 
This  could  be  the  cause  for  some  of  the  error  in  the  calculation  of  the  rms  roughness  by  optical  means.  The  surfaces 
are  close  enough  to  having  a  Gaussian  height  distribution,  however,  that  the  procedure  gives  reasonable  estimations 
of  the  roughness. 

4.  Other  roughness  models 

We  have  now  seen  that  it  is  possible  to  model  certain  rough  surfaces  with  the  one-dimensional  Gaussian  height  FADF 
and  to  obtain  roughness  values  in  this  way,  but  we  have  also  seen  that  this  model  is  often  inadequate  for  real  surfaces. 
We  may  rightfully  ask  if  there  are  other  roughness  models  which  we  may  use. 


4.1.  Isotropic  Gaussian  height  model 


Another  one  of  the  most  common  FADFs  is  the  isotropic  Gaussian  height  model.  This  fimction  describes  surfaces 
which  have  a  Gaussian  distribution  of  height  values  in  all  directions.  This  was  also  proposed  and  studied  by  Beck¬ 
mann  and  Spizzichino;  its  limiting  case  is  die  distribution  function  of  Torrance  and  Sparrow.  It  has  the  form[4] 


5{0^<Py.)  = 


P^ii(cosep^ 


(19) 


Let  us  now  find  the  reflected  radiance  from  a  surface  patch  with  this  FADF.  As  before,  we  may  use  equation  (6)  to 
determine  this  value,  although  we  now  integrate  over  txith  angular  variables.  The  reflected  radiance  will  therefore  be 
given  by 


/’/(cos0^os0j.+  sin0^in0^cos(tp^-(p^))e  ^  sin0^0^qy 
^  ^  fijii {cosQp*cosQ^ 


(20) 


where  0^  and  0^  are  now  both  functions  of  <p^.  The  integral  over  0^  of  the  first  term  in  equation  (20)  is 

-(tanBJ^ 
r-le  ^  . 

I - 5 - ‘“ff 

We  may  integrate  the  second  term  by  parts  over  0^  to  obtain 

J-(tan0^  '  +€//(tanep)d(p^ 

so  that 


't'n 

r  I  —  (tan0..)^  —  (tanOj)^ 

=  j  ^{-e  ^  -I- tan0^cos  (<p^-(pp  (tan0ye  "  -i- e//(tan9p ) ) 


dtp. 
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'/D 
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where  6^  and  9^  are  functions  of  <p^.  We  will  evaluate  this  integral  numerically.  We  are  currently  using  this  equa¬ 
tion  to  model  sandblasted  surfaces,  and  are  testing  the  appropriateness  of  the  model. 

4.2.  General  one-dimensional  FADF 


The  one-dimensional  Gaussian  height  model  and  the  isotropic  Gaussian  height  model  cannot  describe  all  interesting 
surfaces,  despite  the  hopes  of  some  machine  vision  researchers.  Consider,  for  example.  Figure  2S,  which  shows  a  typ- 


column 

Figure  25:  Edge  profile  for  a  ground  surface 

ical  scanline  from  a  metal  plate  which  was  finished  by  grinding.  Notice  that  the  edge  profile  is  almost  linear.  The  facet 
area  distribution  function  for  this  surface  is  obviously  very  different  from  that  proposed  by  Beckmann  and  Spiz- 
zichino,  and  our  current  method  would  therefore  be  of  questionable  value.  It  might  be  possible,  however,  to  obtain 
information  about  the  roughness  of  surfaces  such  as  this  one,  which  are  not  fit  well  by  the  Gaussian  height  models. 

The  method  for  estimating  roughness  which  we  have  analyzed  in  diis  paper  finds  a  single  roughness  parameter  for  die 
surface  from  all  the  points  in  the  edge  profile.  This  is  a  wise  plan  of  attack,  since  it  is  an  overdetermined  system,  and 
is  therefore  fairly  robust  with  respect  to  noise.  We  have  noted,  however,  that  for  one-dimensional  FADF’s  with  certain 
source  and  object  configurations,  the  derivative  of  the  edge  profile  is  closely  linked  to  the  distribution  function. 
Therefore,  rather  than  propose  a  model  of  the  rough  surface  and  find  one  roughness  measure  as  we  do  now,  we  h(^ 
to  be  able  to  obtain  the  shape  of  the  FADF  from  the  derivative  of  the  edge  profile.  Because  there  are  more  degrees  of 
freedom  in  the  output  of  this  method  than  there  are  in  the  output  of  our  present  method,  we  expect  the  derivative 
method  to  be  more  sensitive  to  noise.  The  method  will  also  be  more  sensitive  to  noise  because  the  derivative  operator 
acts  as  a  high  pass  filter.  Nonetheless,  our  initial  experiments  at  obtaining  one-dimensional  FADFs  from  reflected 
edges  are  promising,  and  it  may  soon  develop  into  a  tool  applicable  to  a  general  vision  system.  Figure  26  to  Figure  30 


show  the  FADFs  found  with  this  method  for  the  five  milled  stainless  steel  plates. 
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Figure  26 
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Figure  27 
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FADF  for  7.5  microinches  sample 

Figure  29 
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5.  Implications  for  General  Msion 

The  analysis  presented  above  shows  that  our  method  of  roughness  estimation  poses  no  threat  of  supplanting  current 
methods  of  roughness  measurement.  It  could,  however,  prove  useful  in  general  purpose  vision  systems.  The  ability  to 
analyze  the  reflections  of  step  edges  can  be  \ery  useful  in  artificial  environments,  where  edges  are  ubiquitous,  and 
control  over  the  environment  is  limited.  Additionally,  the  theory  of  edge  detection  and  localization  is  well  developed, 
so  the  initial  steps  in  such  a  vision  system  would  be  based  on  proven  methods. 

Such  a  system  would  begin  by  identifying  the  regions  in  an  image  tqrpropriate  for  this  analysis.  It  would  then  look  for 
the  reflected  images  of  edges  in  each  region.  The  sharpest  reflected  edges  in  the  image  would  probably  come  fitmi 
step  edges  in  the  environment,  and  therefore,  it  would  focus  the  most  attention  on  these  reflected  edges.These  would 
also  have  the  greatest  high  frequency  content,  and  therefore,  diis  search  would  be  facilitated  by  a  frequency  space 
representation  of  the  image,  such  as  the  spectrogram  or  local  Fourier  transform.  Finally,  the  system  would  use  fire 
previously  derived  methods  to  estimate  the  rou^ess  of  the  surface  in  the  given  region. 

In  order  to  construct  a  system  such  as  this,  however,  there  are  several  shortcomings  which  must  be  noted.  The  most 
obvious  is  the  assumption  of  a  known  source  position.  The  position  of  the  source,  i.e.,  the  values  of  9^  and  6^  for 

every  position  on  the  object,  influences  the  reflected  step  edge  profile  drastically.  Hgure  31  shows  die  reflected  edge 
profiles  of  the  same  surface  with  two  different  source  positions.  The  two  curves  appear  quite  different  and  could  be 
mistaken  as  coming  from  surfaces  with  different  roughness  values.  Conversely,  it  is  possible  to  generate  similar 
reflected  edge  profiles  for  surfaces  with  quite  different  roughness  values.  Hgure  32  shows  a  simulation  of  the 
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Figure  32:  IWo  simulated  edge  pn^es.  The  rou^mess 
values  and  source  positions  are  different  for  each. 

positions  for  each  surface,  the  curves  look  very  similar.  In  real  data,  noise  could  easily  erase  any  differences  between 
these  curves.  If  we  consider  equation  (IS),  we  find  part  of  the  problem.  In  many  cases,  the  error  function  terms  in  this 
equation  donunate  the  other  terms,  and  thus,  the  curve  will  be  invariant  for  all  cases  in  which 


/tane-j 'V  /taneUN 

«f/l — J  -  j  equals  a  constant.  Althougli  we  may  not  be  able  to  determine  the  surface  roughness  in 

/tanB^'v  /tan0_'\  tan6^ 

such  cases,  we  may  be  able  to  determine  the  value  of  J  -  erf^ — j  or  — and  — ,  however. 

Another  possible  method  to  determine  roughness  without  a  known  source  position  is  to  estimate  the  source  position  at 
the  same  time  as  estimating  the  roughness.  It  is  also  possible  that  frequency  space  and  wavelet  fcmnulations  might  not 
need  the  knowledge  of  a  source  position. 


Note  that  by  studying  the  derivative  of  the  edge  profile,  we  obtain  information  only  in  the  direction  perpendicular  to 
the  edge,  and  this  limits  the  knowledge  we  can  gain  about  two-dimensional  distribution  functions.  If  there  are  multi¬ 
ple  edges  in  an  image,  however,  we  may  obtain  information  about  the  directional  character  of  the  roughness.  Figure 
33  and  Hgure  34  show  the  scanlines  from  the  same  surface  perpendicular  to  two  different  edges.  The  directional 
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FIgaR  33:  Profile  of  a  step  edge  r^ected  in  a  rough  surface 


character  of  the  toughness  is  obvious  from  these  gnqihs.  By  following  these  lines  of  reasoning,  we  hope  to  improve 
our  methods  of  roughness  estimation  to  make  them  part  of  a  general  machine  vision  system. 


6.  Conclusion 


We  have  shown  that  we  may  estimate  the  roughness  of  metal  surfaces  by  studying  the  reflection  of  step  edges.  The 
method  orders  the  surfaces  the  same  as  would  a  human  observer,  and  this  ordering  appears  related  to  the  rms  rough¬ 
ness  of  the  surfaces.  The  values  obtained  are  similar  to  those  determined  by  other  computer  vision  researchers,  but  are 
very  different  from  those  obtained  by  other  means.  TherefOTe,  the  method  is  probably  impractical  for  precise  mea¬ 
surements  or  quality  control;  it  would,  however,  be  useful  for  a  general  purpose  vision  system.  In  order  to  achieve  this 
goal,  we  must  remove  our  assumption  of  the  known  light  source  position.  We  do  this  in  future  work.  We  also  dis¬ 
cussed  possible  solution  methods  for  surfaces  with  non-Gaussian  height  distributions,  and  the  application  of  such 
methods  to  a  genoal  vision  system. 
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Figure  34:  Profile  of  a  step  edge  rdiccted  in  the  same  rough 
surface,  but  with  the  perpendknlar  orientation 
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9.  Appendix 

We  wish  to  determine  if  the  surfaces  actually  exhibit  a  Gaussian  distribution  of  height  values.  Therefore,  for  each  sur¬ 
face,  we  made  a  histogram  of  the  number  of  occurrences  of  each  height  value  in  the  ten  profiles  for  each  surface. 
These  histograms  appeared  in  Figure  19  to  Figure  25.  As  the  number  of  measured  points  increases,  these  histograms 
should  approximate  a  Gaussian  distribution  mote  and  more  closely.  The  true  form  of  the  histograms  can  be  derived, 
however.  We  first  consider  the  case  of  no  measurement  errors.  Consider  a  discretely  sampled  rough  surface  with  a 
Gaussian  distribution  of  heists.  The  probability  of  a  given  sample  being  in  the  height  range  h  to  h  +  dh  is 

’ 

P  -  L_ - dh 

V2Jia^ 

where  <7^  is  the  standard  deviation  of  the  distribution.  Now,  if  we  have  a  scan  of  the  surface  with  5  sample  points,  the 
probability  that  n  of  them  will  be  in  the  range  from  h  to  h  +  dh  is  given  by  the  binomial  distribution,  or 

probability  =  jp"  ( 1  - 


The  mean  and  variance  of  the  binomial  distribution  are  known  to  be  5P  and  SP(1  -P) ,  respectively. 

If  we  now  allow  the  existence  of  measurement  errors,  we  must  change  the  value  of  P.  The  probability  that  a  mea¬ 
sured  point  has  value  A,  is,  as  before. 


Pi  = 


The  probability  that  the  measuring  instrument  will  record  the  value  of  this  point  as  A  is 


(h-h.y 


dh 


= 
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where  we  have  assumed  that  the  measurement  errors  are  Gaussian  distributed  with  standard  deviation  a  ^  The  prob- 

nt  ^ 

ability  that  the  point  has  value  h^,  but  is  recorded  as  h  is  then  p^q^.  We  sum  these  probabilities  over  all  possible  val¬ 
ues  of  to  determine  the  probability  of  a  point  having  measured  value  in  the  range  from  h  to  h  +  dh.  This 
probability  is 

2of  20^ 

- dh.  ‘—= - 

V2iio.  J2na 

ft  fft 

/  V 

If  «  G^ ,  i.e.  if  the  measurement  errors  are  much  smaller  than  the  standard  deviation  of  the  height  distribution,  as  is 

our  case,  the  second  term  in  the  sum  is  appreciable  only  for  a  small  range  of  values  of  i,  and  the  first  term  is  fairly 
constant  over  this  range.  Therefore,  we  may  approximate  hj  in  the  first  term  by  h,  and  factm*  out  die  first  term  to  find 


^  =  Ip.9,  =  S 


(h-hy 

'  K 


The  sum  is  approximately  equal  to  one,  so  that 

£  * 

/»=£__ - dh 

and  we  may  use  the  results  we  derived  for  the  case  with  no  measurement  errws. 
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